Inverting Time-Dependent Harmonic Oscillator Potential by a Unitary 
Transformation and a New Class of Exactly Solvable Oscillators 



Ali Mostafazadeh* 



Theoretical Physics Institute, University of Alberta, 



^ , Edmonton, Alberta, Canada T6G 2J1. 

Q\ _ 

_^ 

November 1996 

> 
O 

f***- ' Abstract 

A time-dependent unitary (canonical) transformation is found which maps the Hamiltonian for a harmonic 
oscillator with time-dependent real mass and real frequency to that of a generalized harmonic oscillator with 
time-dependent real mass and imaginary frequency. The latter may be reduced to an ordinary harmonic oscillator 



> 



' by means of another unitary (canonical) transformation. A simple analysis of the resulting system leads to the 

identification of a previously unknown class of exactly solvable time-dependent oscillators. Furthermore, it is 
shown how one can apply these results to establish a canonical equivalence between some real and imaginary 
frequency oscillators. In particular it is shown that a harmonic oscillator whose frequency is constant and whose 
mass grows linearly in time is canonically equivalent with an oscillator whose frequency changes from being real 
■ to imaginary and vice versa repeatedly. 

The solution of the Schrodinger equation, Hip = itp, for a harmonic oscillator with time-dependent mass m and 
frequency u>, i.e., 

£3 ! 1 m{t)uj 2 (t) o 

has been the subject of continuous investigation since late 1940's, jl], 0,0,0). The main reason for the interest 
in this problem is its wide range of application in the description of physical systems. Although by now there 
exist dozens of articles on the subject, a closed analytic expression for the time-evolution operator is still missing. 
Recently, Ji and Kim [Q showed that using the Lewis-Riesenfeld method Q one can construct an invariant operator 
in terms of the (two independent) solutions of the classical dynamical equations: 

i[m[t) j t x c (t)} + m(t)uj 2 {t)x c (t) = , (2) 

and therefore reduce the solution of the Schrodinger equation to that of Eq. (|J), for u)(t),m(t) 6 IR. The case 
where the frequency ui is imaginary has been considered only in the time- independent case g| . 

The purpose of this note is to study the implications of the recently developed method of adiabatic unitary 
transformation of the Hilbert space [0 for this problem. The basic idea of this method is to use the inverse of 
the adiabatically approximate time-evolution operator to transform to a moving frame. This transformation has 
proven to lead to some interesting results for the system consisting of a magnetic dipole in a changing magnetic 
field. The analogy between the dipole system and the time-dependent harmonic oscillator is best described in terms 
of the relation between their dynamical groups, namely SU{2) and SU(1, 1), 0. 
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Let us first concentrate on the real frequency case, u>(t) € K. Then the Hamiltonian ([j]) and its eigenvectors 
\n:t) can be expressed in terms of the creation (a') and annihilation (a) operators, namely 

H(t) = Lo(t)[a)(t)a{t) + 1/2] , (3) 
'nit) = -±=a^(t)\0;t) , (4) 



i n 
where 

a(t) := ±= (e^x + ie~<%) , «(t) := ± ln[m(t)w(t)] , (5) 

and |0;t) corresponds to the ground state at time t, i.e., a(t)\0;t) — 0. Similarly to the time-independent case the 
eigenvalue equation 

H(t)\n;t) =E n (t)\n;t) , (6) 

is satisfied for E n (t) — E n [w(t)] = <*>(£) (n+ 1/2). Therefore the eigenvalues are non-degenerate and the adiabatically 
approximate time-evolution operator may be expressed as: 

U (t):=J2e ia ^\n;t)(n;0\, (7) 

n 

where 

<*»(*) := «„(*)+ 7n(*) , *n(*) := - / E n {t')dt' 7n (t) := t / A nn (t')^', (8) 

Jo Jo 

A ro „(i) := (m;t|-||n;t) . (9) 

Following the ideas developed in Ref. let us next use Uq 1 ^) = U^(t) to perform a unitary transformation 
of the Hilbert space. Recall that under a general unitary transformation \ip(t)) — > = W(t)|^>(i)), the 

Hamiltonian H and the corresponding time-evolution operator transform according to 

H{t)^H'(t) = U{t)H{t)U\t)-iU{t)U\t) , (10) 
U(t)->U'(t) = U{t)U{t)U\0) . (11) 



Substituting Uq for in (10), one can show, ||, that 

= "< E e- i[Qm(t )- a " W U ro „(i)|m; 0)(n; 0| . (12) 

In order to express the transformed Hamiltonian if' in a closed form, one must first compute A mn . This can 
be done directly by substituting Eq. (^) in (|]) or indirectly using the identity 

(m-t\ dH ^\n ■ t) 

which is obtained by differentiating both sides of Eq. (^|) and making use of the orthonormality of \n\ t). The latter 
method turns out to be much simpler. It yields 

Amn(t) = —-- V n(n - l)6 m , n -2 - V m(m - l)8 m -2,n • (13) 
Here use is made of the well-known relations 

a(t)\n;t) =y/n\n-l;t) , a^(t)\n;t) = Vn+ 1 \n + l;t) , and a(t) = k(t)at(t) . (14) 
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Note that Eq. ( p"3| ) is also valid for m — n since due to the fact that |n; t) can be chosen to be real, A nn (t) = 0. In 
particular 7„(t) = 0. Hence, a n (t) — 5 n (t) = (n + l/2)S(t), where 

S(t) :=- [ uj{t)<1t . (15) 
Jo 

Substituting the expressions for a n and A mn in Eq. ([l2]) and performing the summation over m, one finds 
-ik(t) 



H'(t) = 



-ik(t) 



2 

«(*) 
2 



5^ [Vn(n - l)e 2i5 W|n - 2; 0){n; 0| - y/(n + l)(n + 2) e - 2i4(t) |ra + 2; 0)(n; 0| 

n 



^(0) 



-2i5(i)„f2 



a^(0) 



{sin[2<5(0](e- 2K( V - e 2K(0) x 2 ) + cos[2<5(i)](xp + pz)} 



(16) 



Next consider the Schrodinger equation in the transformed frame: H'(t)\ip'(t)} = i\ip'(t)). The presence of a total 
derivative k on the right hand side of Eq. ([if]), suggests a redefinition of the time t — > t' := «(t). Note that for 
k{t) = 0, i.e., Lu(t) = c/m(t) for some constant c G IR, H'{t) vanishes identically and the adiabatic approximation 
is exact. Furthermore, for k < 0, one can consider the time-reversed system for which k > 0. The time-evolution 
operator for the original system is obtained from that of the time-reversed system by inversion. Therefore, without 
loss of generality, one can assume k(t) > 0. The latter allows for the above-mentioned redefinition of time t — * t' . 
This leads to the Schrodinger equation h(t')\ip' (t')) = i\ip'(t')) for the Hamiltonian 



h(t') := ^{sm[25{t')}{e- 2Ko p 2 -e 2Ko x 2 )+cos[2S(t')}{x P +px)} , 



1 2 

2m'{t') P + 



m '{t')J 2 {t') 



1 



x 2 + -0(t')(xp + px) , 



where 



5(f) 



S(t(t')) , k := K (0) , m'(i') 



sin[2£(t')] 



J{t') := ism[25{t')}, 0(1?) := cos[2<5(t')] = ^l+u' 2 (t'). 
Hence, the transformed Hamiltonian is a generalized harmonic oscillator JtJ: 

MO 

a(t') 



1 



a(t')p 2 + /3(t')(xp + px) + -i(t')x 2 } 
1 (t') := m'(t')J 2 (t') 



2 

l/m'(t') 



(17) 

(18) 
(19) 

(20) 



with a real mass m'(t') and an imaginary frequency w'(t') = ie 2Ka /m'(t'). Note that as a result of the adiabatic 
unitary transformation and redefinition of time, the two arbitrary functions m(t) and oj(t) have been reduced to a 
single function namely d(t'). 

It is well-known that one can transform the generalized harmonic oscillator to an ordinary harmonic oscillator 
by the time-dependent canonical transformation, |Q , 

"/?(*')" 



x — > x 



p^p + 



a(t') 



This leads to the Hamiltonian 



h>(t>) 



a(t')p 2 



l(t') 



f3 2 (t') d (0(f) 
a(t') dt \a(t') 



1 2 , m'(t')ri*(t') 2 
2m'(t') P + 2 X 



where 



Q'(t') :=W-1 



i[2<J(t')] tft 7 



(21) 



(22) 



(23) 
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Note that W can be real or imaginary depending on the form of 5(t'). Clearly for il' = 0, the problem reduces to 
a free particle with a variable mass whose solution can be exactly given. In terms of the original functions, the 
condition f2' = is expressed as 



m(t) 



m 
w(t) 



tan 



Lo(r)dT 



(24) 



where mo is a real constant. 

Eq. ( p4| ) determines a new class of exactly solvable cases which is the analog of the exactly solvable magnetic 
dipole Hamiltonians obtained in Q. The only difference is that here an additional canonical transformation ( pl| ) 
is also performed. It is not difficult to see that this transformation corresponds to the action of the unitary 
operator U'(t') := exp[i/3(t')x 2 / (2a(t))] on the Hilbert space. In fact, the Hamiltonian h! can be obtained from 
h by substituting h for H and U' for U in Eq. (|ic|). Similarly in view of Eq. ([Tl|), one has the following relation 
between the evolution operators 

u(0 = e~ l — x u'{t') e*^> x , (25) 



where it(t') and u'{t') are the evolution operators for the Hamiltonians h(t') and h('(t'), respectively. A further 
application of Eq. (nl]) leads to the expression 



U(t) = U (t)u(t'(t)) , 



(26) 



for the time-evolution operator associated with the original Hamiltonian (Q). Here Uo(t) is the adiabatically 
approximate expression for the time-evolution operator given by Eq. (^). 

For the cases where the condition (Q) holds, u'(t') is the evolution operator for a free particle with a variable 
mass m'(t'), i.e., 

' -i /"*' dr 

If f2' does not vanish but satisfies £l'(t') — £l' /m'(t') for some real constant fl' , then the Schrodinger equation 
for h! is still exactly solvable. In fact in this case 



u'{t') — exp 



(27) 



u'(i') = exp 



dr 
m'(r) 



(p 2 + n 2 x 2 ) 



The condition Q,'(t') — Q,' /m'(t') is a generalization of Q'(t') = 0. It corresponds to a larger class of exactly solvable 
time-dependent harmonic oscillators. In terms of the original parameters m and uj this condition is expressed as: 



m(t) 



(28) 



where mo and C := 1 



/fig are positive constants and 



fit) ■-- 



9{t) 



z{t) := 1 



C-i 



Lj{r)dT 



Eq. ( p4| ) is a special case of ( pq ) where — > 0. 

Finally I would like to emphasize the following points: 

1) For the case that fl' is real, one can repeat the above analysis by replacing H of Eq. (|l|) by b! of Eq. (p2|). In 
principle this may lead to yet other exactly solvable cases. If the iteration of this procedure yields oscillators 
with real frequency at each step, then it can be repeated indefinitely. This leads to a product expansion for 
the time-evolution operator which is analogous to what is called an adiabatic product expansion in Ref. fl. 
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2) The above analysis also indicates that the time-dependent harmonic oscillators whose mass M and frequency 
O are related according to 

/ m" 

n=i/-l ; (29) 

play a universal role. This is because as shown above the problem for the most general real frequency harmonic 
oscillator (Q) can be reduced to this case by means of a series of unitary (canonical) transformations. Eq. ( p9|) 
is obtained from (|23|) by expressing the right hand side of ( ffit ) in terms of the mass. The parameter Mo is 
an arbitrary positive constant corresponding to e . 

3) It is tempting to seek applications of the known results for the real frequency oscillators to the time-dependent 



imaginary frequency oscillators satisfying (29). In view of the above analysis, there is a class of imaginary 

frequency oscillators of this form which are canonically equivalent to some real frequency oscillators. One 

must however be aware that by performing the canonical transformation described in this article in the reverse 

order, one might not be able to transform an arbitrary imaginary frequency oscillator to a real frequency 

one, even if it satisfies Eq. ( ^9| ) . In this case the transformation associated with the adiabatic approximation 

would not be the same as the one obtained above. Nevertheless, the above scheme is consistent in the 

sense that if the frequency of the oscillator obtained by transforming back an imaginary frequency oscillator 

satisfying ( p9| ) turns out to be real then one has the desired result. Otherwise, the method fails to transform 

the imaginary frequency oscillator to a real frequency one. A simple example of a case where Eq. ( |2^ ) is 

satisfied but the transformation to a real frequency oscillator is not possible is the time-independent oscillator 

with M = const, and f2 = i. Another example is the case where the mass is decaying in time £' according 

to M = m'(t') = M (l + e^ 4 ') with fj, > 1 and the frequency Q = £}'(£') is given by Eq. i.e., 

O 2 = — 1 + /i(l + 2e Mt )~ 1 / 2 . This system is particularly interesting since the potential V := MQ 2 /2 changes 
2 l 

sign at t' = k* := ^ ln( 7 ~ ). It is positive for £' < and negative for t' > k*. Hence, at t' = ft* the energy 
spectrum undergoes a 'phase transition.' However, it is not difficult to see that transforming this system to 
a real frequency oscillator is not possible. This is because making the necessary canonical transformations, 
one obtains an oscillator with the mass m(£) and frequency u(t) satisfying: 

[m(tMt)}- 2 /» = 1 - 1 (30) 
sin 2 J Q Lu(T)dT 

Thus one should seek functions m : [0, oo) — > H + and uj : [0, oo) — > K + which satisfy Eq. ( ^0| ) and are 
positive for £ € [0,T] for some T > £*, where f» is defined according to k(£») = ft*. That such functions 
do not exist can be directly inferred from the fact that for £ — > the right hand side of Eq. ( |30| ) tends to 
— oo, whereas the left hand side remains positive. A simple example of a real frequency oscillator which is 
canonically equivalent to an oscillator whose frequency fluctuates between real and imaginary values is 

m(£)=mo + /x£, w(t)=uo, (31) 

where mo, /i, and loq are positive constants. Applying the canonical transformations introduced in this article 
to this oscillator, one arrives at another canonically equivalent oscillator with mass m'(t') and frequency Q'(t') 
given by 

m'(t')= , m °"° 71 , tt' 2 (£') = -l ~t j (32) 

sin -(m o;o-e 2t ) /xsin -(m uj - e 2t ) 

where £' = k(£) = ln[wo(wo + fit)]/2, and use is made of Eqs. ( ^8| ) and (]23|). Fig. 1 shows a plot of £1 2 as 
a function of £', for mo = (J, = ujq = 1. As seen from this plot, fl 2 changes sign repeatedly. However, the 
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Figure 1: Plot of 2 as a function of t' for the oscillator of Eq. ( |32| ) with mo = // = u>o = 1. 

original oscillator ( pl| ) is a very simple time-dependent oscillator with positive real mass and frequency. In 
fact, one might try to apply the results of Refs. ||, Q to obtain exact solution of the Schrodinger equation 
for this oscillator. This would immediately lead to the exact solution of the Schrodinger equation for the 
oscillator (p3), at least for the periods of time during which m! and il' are continuous functions of t' . This is 
particularly interesting since as fl 2 changes sign from positive to negative, the spectrum of the corresponding 
oscillator ( |32| ) changes from being discrete to continuous and vice versa, while the spectrum of the canonically 
equivalent oscillator (EHl) remains always discrete. 
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